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Abstract 

Recent Wilkinson Microwave Anisotropy Probe (WMAP) data confirm 
the Cosmic Microwave Background (CMB) quadrupole anomaly. We further 
elaborate our previous proposal that the quadrupole power can be naturally 
suppressed in axis-symmetric universes. In particular, we discuss in greater 
detail the CMB quadrupole anisotropy and considerably improve our anal- 
ysis. As a result, we obtain tighter constraints on the direction of the axis 
of symmetry as well as on the eccentricity at decoupling. We find that the 
quadrupole amplitude can be brought in accordance with observations with 
an eccentricity at decoupling of about 0.64 x 1CP 2 . Moreover, our determi- 
nation of the direction of the symmetry axis is in reasonable agreement with 
recent statistical analyses of cleaned CMB temperature fluctuation maps ob- 
tained by means of improved internal linear combination methods as Galactic 
foreground subtraction technique. 
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1 Introduction 

The Cosmic Microwave Background (CMB) has had a profound impact on mod- 
ern cosmology and greatly improved our understanding of the universe. The CMB 
angular power spectrum is indeed very sensitive to the origin and evolution of the 
cosmic density fluctuations. 

The temperature fluctuations of CMB are observed at the level of AT/ (T) ~ 
10~ 5 [lj. The high resolution data provided by the Wilkinson Microwave Anisotropy 
Probe (WMAP) [21 El HI [5] confirmed that the CMB anisotropy data are in striking 
agreement with the predictions of the simplest inflation model. 

However, the 3-years WMAP data (WMAP3) display at large angular scales 
some anomalous features. The most important discrepancy resides in the low 
quadrupole moment, which signals an important suppression of power at large 
scales, although the probability of quadrupole being low is not statistically com- 
pelling. The problem consists into the fact that the power of the quadrupole is 
substantially reduced with respect to the value of the best-fit A-dominated cold 
dark matter (ACDM) standard model. If this discrepancy turns out to have a cos- 
mological origin, then it could have far reaching consequences for our understanding 
of the universe and in particular for the standard inflationary picture. This peculiar- 
ity emerged since 1992 when the first data of the differential microwave radiometer 
(COBE/DMR) appeared [TJ. Since then, in 2003 (WMAP) and in 2006 (WMAP3) 
this behavior was confirmed. In fact, the WMAP3 data give a quadrupole power of 
211 /iK 2 , while the expected value in A-dominated cold dark matter model is about 
1252 /j,K 2 . 

In the last years, the "smallness" of CMB quadrupole has been subject to very 
intensive studies because it may signal a non-trivial topology of the large scale ge- 
ometry of the universe [6j. Indeed, several possibilities have been advanced in the 
recent literature to understand the suppression of the quadrupole power [7J E] (for 
other large scale anomalies in the angular distribution of CMB see Ref. [H]). 

Recently [10] we showed that, allowing the large-scale spatial geometry of our 
universe to be plane-symmetric with eccentricity at decoupling of order 10 -2 , the 
quadrupole amplitude could be drastically reduced without affecting higher multi- 
poles of the angular power spectrum of the temperature anisotropy. Remarkably, 
the "ellipsoidal" (or "eccentric" ) universe has been considered also in non-standard 
cosmological models, such as braneworld cosmology [TT] . 

There are different mechanisms which could induce a planar symmetry in the spatial 
geometry of the universe [12] • Among these, the most interesting examples include 
a cosmic domain wall, a cosmic string and an almost uniform cosmic magnetic field. 
In particular, the cosmic magnetic field seems to be of relevance since several ob- 
servations suggest the presence of magnetic fields correlated on cosmic scales [13] 
(for a recent study of the effects of cosmic magnetic fields on large-scale structures 
of the universe see Ref. [13], while for influence of cosmic fields on the expansion 
of the universe see Ref. [IS]). There are several possible mechanisms to generate 
cosmological magnetic fields. For instance, magnetic fields might be produced dur- 
ing [TBI [TT] or after [18] the inflation era. 
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In this paper, we further discuss the proposal advanced in Ref. [10] where we 
showed that cosmic microwave background quadrupole power can be naturally sup- 
pressed in plane-symmetric universes. In particular, we discuss in greater detail the 
CMB quadrupole anisotropy and considerably improve the analysis presented in our 
previous paper [10J. As a result, we obtain tighter constraints on the direction of 
the axis of symmetry as well as on the eccentricity at decoupling. 

The plan of the paper is as follows. In Sec. 2 we discuss the Einstein's equations 
for cosmological models with planar geometry, and we describe the mechanisms to 
generate the eccentricity in the universe expansion by means of magnetic fields, do- 
main walls or cosmic strings; Sec. 3 deals with the analysis of CMB anisotropies 
including the asymmetric contributions to the temperature anisotropy. In Sec. 4 
we discuss the constraints on cosmic magnetic fields arising from primordial nucle- 
osynthesis the large scale structure formation. Finally, we draw our conclusions in 
Sec. 5. Some technical details are relegated in the Appendix. 



2 Cosmological models with planar symmetry 

We are interested in cosmological models with planar symmetry. The most general 
plane-symmetric line element [19] is: 



ds 2 = dt 2 - a 2 {t){dx 2 + dy 2 ) - b 2 {t) dz 1 



(2.1) 



where a and 6 are the scale factors. The metric 112. ip corresponds to considering the 
xy-plane as a symmetry plane. 

The non-zero Christoffel symbols corresponding to the metric (12.11) are: 



1 li 



1 22 



aa, r 33 



66, rl, 



r 2 

1 02 



a/ a, r 3 03 = b/b, 



[2.2) 



where a dot indicates the derivative with respect to the cosmic time. The non-zero 
Ricci tensor components turn out to be: 



Rn 



R\ 




Rr, = 



a b 
a b 



o ao 



(2.3) 
(2.4) 
(2.5) 



The most general energy-momentum tensor consistent with planar symmetry is 

Tt = diag(p, -p\\,-p\\,-p±). (2.6) 
It can be made up of two different components: an anisotropic contribution, 

{TaT v = diag (p A , -pf , -pjf , -pi), (2.7) 
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which induces the planar symmetry -as, for example, a uniform magnetic field, a 
domain wall, or a cosmic string-, and an isotropic contribution, 



{TiY u = diag (p 1 , -p 1 , -p 1 , -p 1 ), (2.8) 

such as vacuum energy, radiation, matter, or cosmological constant. Exact solutions 
of Einstein's equations for different kind of plane-symmetric plus isotropic compo- 
nents can be found in Ref. [12J. 

Taking into account the above energy-momentum tensors, the Einstein's equations 

R»u - \g^R = 8vrGT^, (2.9) 

read 

(-Y + 2*£ = 87rG(/ + /), (2-10) 

a h ah „. T 4x , 

- + u + ~u = ~^G(p ! +pf , 2.11 
a b a b 11 

2-+(-) =-8nG(p I + pf). (2.12) 
a \a J 

In the following we shall restrict our analysis to the case of matter-dominated uni- 
verse {p 1 = 0) filled with an anisotropic component given by a uniform magnetic 
field (directed along the z- axis), or a cosmic domain wall (whose plane of symmetry 
is the xy-plane), or a cosmic string (directed along the z-axis). 
Magnetic fields have been observed on a wide range of scales. In particular, they 
have been detected in galaxies, galaxy clusters, and also in extra-galactic structures 
(for recent reviews on cosmic magnetic fields, see Ref. p3]). It is reasonable to as- 
sume that the actual observed magnetic fields have a primordial origin. We assume 
that a (almost) uniform magnetic field pervades our universe, though it is not ex- 
cluded that such a field may have a more complicated structure on small scales. 
Further examples of anisotropic components are given by cosmic topological de- 
fects [2D]. As it is well known, phase transitions in the early universe can generate 
domain walls (cosmic strings), which are two-dimensional (one-dimensional) defects 
that originate when a discrete (axial or cylindrical) symmetry is broken. 

For the three cases discussed above, the energy-momentum tensor take on the 
form 

(T B )^ = p fl diag(l, -1,-1,1), (2.13) 

(T w y v = p w diag(l, 1,1,0), (2.14) 

(T s )^ = p s diag(l, 0,0,1), (2.15) 

where p#, p w , and p s , are the magnetic, wall, and string energy density, respectively. 
Moreover, we assume that the interaction of anisotropic components with matter 
is negligible (in the case of magnetic fields, this corresponds to taking into account 
that the conductivity of the primordial plasma is very high [1 3J ) . In this case, the 
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anisotropic component of the energy-momentum tensor is conserved, (T^^ = 0, 
so that we have 

p A + 2-(p A +pi) + Up A +pi) = 0. (2.16) 
a " b 

Let us introduce the eccentricity 



6V « -\A-(^ ^ 

and normalize the scale factors such that a (to) = b(to) = 1 at the present time to- 
The first definition of eccentricity applies to the cases of uniform magnetic field or 
string, while the second one to the case of domain wall. 

In this paper, we restrict our analysis to the case of small eccentricities (that is 
we consider the metric anisotropies as perturbations over the isotropic Friedmann- 
Robertson- Walker background). In this limit, from Eqs. (I2.10I) - (I2.12I) . we get the 
following evolution equation for the eccentricity: 

^^ + 3H{ee) = ±87iG{pf-pi), (2.18) 

where the plus sign refers to the cases of magnetic field and string, while the minus 
sign to the case of domain wall. Here H = a /a is the usual Hubble expansion 
parameter for the isotropic universe. In the matter-dominated era, it results a(t) oc 
t 2 / 3 , so that H = 2/(3t). 

Moreover, to the zero-order in the eccentricity, from Eq. (12.161) it follows that the 
energy densities (and pressure) scale in time as ps oc a~ 4 , p w oc a" 1 , and p s oc a -2 , 
for the three cases, respectively. 

The solution of Eq. ( 12. 18ft for the magnetic case is 



e 



2 _ an(°) 



8fig ) (l-3a- 1 + 2a- 3/2 ), (2.19) 



where = Pb(£q)/Pc?, and p£r = 3H^/8ttG is the actual critical energy density. 
At the decoupling, t = t dec , we have e\ ec ~ 16 where ed ec = e(tdec) and 

Zdec — 1088 is the red-shift at decoupling [4J. Accordingly, we get: 



,_2 i4 B 



1/2 



or 



10 "Mipl . ( 2 - 2 0) 



lO' 2 ^' 1 — (2 21) 

U 10- 8 G' 1 J 



where Bq = B(to) and h ~ 0.72 0] is the little- h constant. 

If, for instance, we assume for the present cosmological magnetic field strength the 
estimate B ~ 5 x 10~ 9 G, which is compatible with the constraints analyzed in 
Ref. [2~TI . we get an eccentricity at decoupling of order e& ec ~ 10~ 2 . 



5 



In the cases of domain wall and cosmic string, integrating Eq. (I2.18p . we find 

e 2 = ^W(3a 2 + 4a- 3 / 2 -7) (2.22) 



and 



e 2 = -nf(3a + 2a- 3 / 2 -5), (2.23) 
5 

respectively, where £1$ and are the actual energy densities, in units of pc?, 
associated to domain wall and cosmic string. 

From the above equations, we can estimate the eccentricity at decoupling in terms 
of wall and string energy densities at present time: 

1/2 

(2.24) 

and 

/ o(°) V /2 

(2.25) 

Vice versa, since the analysis of the CMB radiation constrains the value of the 
eccentricity at decoupling to be less than 10~ 2 (see Ref. [10J and next Section), one 
gets an upper limit on the values of the energy density of cosmic defects stretching 
our universe of order < 1CT 7 . 

3 CMB Quadrupole Anisotropy 

Let us begin by briefly discussing the standard analysis of the CMB temperature 
anisotropics [22]. First, the temperature anisotropy is expanded in terms of spherical 
harmonics: 

AT(M) 




(T) 

> ' 1=1 m=-l 

After that, one introduces the power spectrum: 



a l™Ylm(0A)- (3-1) 



AT Z _ 1 /(/ + 1) ^, 

(T) "y27r 2/ + 1 ^ lml ' 1 } 

that fully characterizes the properties of the CMB anisotropy. In particular, the 
quadrupole anisotropy refers to the multipole £ = 2: 

S-^, (3.3) 

where (T) ~ 2.73K is the actual (average) temperature of the CMB radiation. The 
quadrupole problem resides in the fact that the observed quadrupole anisotropy is 
in the range (see Table 1): 

(AT 2 ) 2 bs ^(210- 276)/iK 2 , (3.4) 
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while the expected quadrupole anisotropy according the ACDM standard model is: 



(AT 2 )j 2 ~ 1252/iK 2 . (3.5) 

If we admit that the large scale spatial geometry of our universe is plane-symmetric 
with a small eccentricity, then we have that the observed CMB anisotropy map is a 
linear superposition of two contributions [TU1 123] : 

AT = AT A + ATi , (3.6) 

where ATa represents the temperature fluctuations due to the anisotropic space-time 
background, while ATI is the standard isotropic fluctuation caused by the inflation- 
produced gravitational potential at the last scattering surface. As a consequence, 
we may write: 

O'lm = dim + a lm- (3-7) 

We are interested in the distortion of the CMB radiation in a universe with pla- 
nar symmetry described by the metric (|2.1|) . As before, we will work in the small 
eccentricity approximation. From the null geodesic equation, we get that a pho- 
ton emitted at the last scattering surface having energy E dec reaches the observer 
with an energy equal to E (n) = (E )(l - e\ cc nl/2), where (E ) = E dec /(1 + z dcc ), 
and n = (ni, n-2, n^) are the direction cosines of the null geodesic in the symmetric 
(Robertson- Walker) metric. 

It is worth mentioning that the above result applies to the case of the axis of symme- 
try directed along the z-axis. We may, however, easily generalize this result to the 
case where the symmetry axis is directed along an arbitrary direction in a coordinate 
system (x g , y g , z g ) in which the x g ?/ g -plane is, indeed, the galactic plane. To this end, 
we perform a rotation H = Tl x (&)H z (p + tt/2) of the coordinate system (x,y,z), 
where H z (p + tt/2) and H x (&) are rotations of angles if + n/2 and $ about the z- 
and a;-axis, respectively. In the new coordinate system the magnetic field is directed 
along the direction defined by the polar angles (p). Therefore, the temperature 
anisotropy in this new reference system is: 

AT A _ E (n A ) - (E ) _ 1 

ITT = (Eo) " ~ 2 ec ' w 

where = (Jin) 3 is equal to 

n\(9, 4>) = cos 6 cos 1? — sin 6 sin-i? cos(0 — <p) . (3.9) 

Alternatively, when the eccentricity is small, Eq. (12. ip may be written in a more 
standard form: 

ds 2 = dt 2 - a 2 (t)(5 lj + hij) dx l dx\ (3.10) 
where is the metric perturbation which takes on the form: 

= -e 2 5 i3 5 j3 . (3.11) 
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The null geodesic equation in a perturbed Friedmann- Robert son- Walker metric gives 
the temperature anisotropy (Sachs- Wolfe effect): 

AT 1 f to dh- 

1 dt^nW, (3.12) 



where n l are the directional cosines. Using e(to) = 0, from Eqs. ( 13. lip and ( I3.12j) 
one gets: 

"^y = — 2 e dec n L (3.13) 
which indeed agrees with our previous result. 

It is easy to see from Eq. (13.80 that only the quadrupole terms (i = 2) are different 
from zero: 



7T 



6^ 



2i. - -w=[l + 3cos(2t?)]eL, 



4i = = ^ e~* sin(2tf) eL , 

4 = «- 2 r = -V^ e " 2i " sin¥e ^- (3 - 14) 

Consequently, the quadrupole anisotropy is: 

Q A = -^=e 2 dec . (3.15) 

Since the temperature anisotropy is a real function, we have a^_ m = (— l) m (a/ jm )*. 
Observing that a£_ m = (-l) m (<J* [see Eq. jSUp], we get a}'_ m = (-l) ro (a£j*. 
Moreover, because the standard inflation-produced temperature fluctuations are sta- 
tistically isotropic, we will make the reasonable assumption that the a\ m coefficients 
are equals up to a phase factor. Therefore, we can write: 

4i=-(4-i)*=yf e ^ 2 Qi> ( 3 - 16 ) 

4 2 = (4,- 2 )* = yf e^d, 

where < 0, < 27r are unknown phases, and 

Qi ~ 13 x 1(T 6 . (3.17) 

Taking into account Eqs. (13.141) and (I3.16p . and Eqs. (13. ip . (13.21) and (13.71) . we get 
for the total quadrupole: 



Q 2 = Ql + Qi - 2fQ A Q h (3.18) 
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Table I: The cleaned maps SILC400, WILC3YR, and TCM3YR. Note that the 
values of a 2m in this table correspond to the values of a 2m given in Refs. [2U El ES] 
divided by (T) ~ 2.73K. The values of the angles 02, and 3 are in degrees. 



MAP 


m 


Re[a 2m ]/l(r 6 Im 




/io- 6 


02 03 


(AT 2 )7/iK 2 







2.75 


0.00 






SILC400 


1 


-0.56 


1.77 


107.6 44.2 


275.8 




2 


-6.79 


-6.60 











4.21 


0.00 






WILC3YR 


1 


-0.02 


1.78 




90.6 52.5 


248.8 




2 


-5.28 


-6.89 











1.22 


0.00 






TCM3YR 


1 


0.10 


1.79 


86.8 49.2 


209.5 




2 


-5.45 


-6.32 







where 

/(#, (p ; 0i, 02) 03) = — 7= {2^ [sin?9cos(2(y9 + 3 ) — 2 cosi9cos(y? + 2 )] sin?? 

4y5 

+ [1 + 3cos(2^)]cos0i}. (3.19) 

Looking at Eq. (13.181) we see that, if the space-time background is not isotropic, 
the quadrupole anisotropy can become smaller than the one expected in the stan- 
dard picture of the ACDM (isotropic-) cosmological model of temperature fluctua- 
tions. We may fix the direction of the magnetic field and the eccentricity by solving 
Eq. (13.71) . which is (for I = 2) a system of 5 equations containing 5 unknown pa- 
rameters: edec, V 9 , 02, and 3 . Note that it is always possible to choose a 20 real, 
and then 0i = 0. 

To solve Eq. (13. 7\i for £ = 2, we need the observed values of the a 2m 's. We 
use the recent cleaned CMB temperature fluctuation map of the WMAP3 data 
obtained by using an improved internal linear combination method as Galactic fore- 
ground subtraction technique. In particular, we adopt the three maps SILC400 [24J, 
WILC3YR [5], and TCM3YR [25]. For completeness, we report in Table 1 the values 
of a 2m corresponding to these maps. We also indicate the values of 2 and 3 [found 
by solving Eq. (I3.16P ] in the case of isotropic universe (a 2m = 0). Moreover, for sake 
of definiteness, we assume that the planar symmetry is induced by a cosmological 
magnetic field with strength B at the present time. 

Numerical solutions of Eq. (13. 7p . referring to the three maps, are given in Tables 
2, 3, and 4, respectively. In the Appendix, we will show that the system (13.7P admits 
at most 8 independent solutions. However, due to the particular values taken by the 
a 2m for the three different maps, we find that to the map SILC400 it corresponds 8 
independent solutions, while to the maps WILC3YR and TCM3YR there correspond 
only 4 independent solutions. 

Moreover we observe that, for each independent solution (edec, V 9 , 02, 03) shown in 
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Table 2: Numerical solutions of Eq. (13.71) obtained by using the map SILC400; the 
values of the angles <p, 2 , and 3 are in degrees. 



edec/10- 2 


■& 


<P 


02 


03 


B /10- 9 G 


0.69 


38.6 


80.1 


100.8 


34.8 


5.0 


0.69 


38.4 


99.2 


83.7 


20.9 


4.9 


0.68 


38.3 


37.2 


139.1 


66.8 


4.9 


0.67 


37.4 


139.3 


47.3 


12.8 


4.8 


0.63 


35.3 


308.9 


43.2 


17.3 


4.5 


0.62 


34.5 


281.2 


73.6 


25.8 


4.5 


0.62 


34.3 


239.8 


122.5 


48.9 


4.4 


0.62 


34.3 


255.2 


104.2 


39.9 


4.4 



Table 3: Numerical solutions of Eq. ( 13. 7ft obtained by using the map WILC3YR; 
the values of the angles <p, 2 , and 3 are in degrees. 



edec/10- 2 







02 


03 


B /10- 9 G 


0.69 


40.6 


91.7 


88.6 


28.0 


4.9 


0.69 


40.6 


88.3 


91.6 


30.8 


4.9 


0.67 


40.0 


139.4 


45.6 


8.6 


4.8 


0.67 


40.0 


40.4 


134.6 


71.0 


4.8 



Table 4: Numerical solutions of Eq. (13.71) obtained by using the map TCM3YR; the 
values of the angles <p, 2 , and 3 are in degrees. 



edec/10- 2 





f 


02 


03 


B /10- 9 G 


0.70 


36.5 


96.6 


83.8 


24.9 


5.0 


0.70 


36.5 


83.2 


95.7 


35.1 


5.0 


0.69 


36.1 


129.4 


54.3 


9.4 


5.0 


0.69 


36.0 


51.6 


124.1 


59.8 


5.0 
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the tables, there exists another one given by (edec,^ — V 9 =t 71 1 02,03), where we 
take the plus sign if ip < n and the minus sign if <p > ir. 

Looking at Table 1 and Eq. (13.1 7p , we see that the values of coefficients | a 2 i | are 
about one order of magnitude smaller then Q\. Assuming Q\ ^> |a 2 i|, we will show 
in the Appendix that an approximate solution for ed ec and d is: 



^dec 



10 Qi 



(3.20) 



and 



d ~ - arccos 
2 



(5VE-VEir)Qi - 5V3 



«20 



3V67rQi 



(3.21) 



From Eq. fl3T20D we also have that e dec ^ 0.64 x 10~ 2 , B ~ 4.6 x 10~ 9 G, and 
Q A ~ (4/y / 30)Qi — 0.7Qi. From Eq. (13. 21 p we get, for the three different maps, 
$ ~ 34°, 36°, 31°, respectively. As one can easily check, the approximate solutions 
Eqs. (I3.20p and (13.211) are quite close to the numerical values. 



50 

b 

40 

30 



WILC3YR 

20 ^ ■ TCM3YR 

= • SILC 400 
10 — 

o Lj — I — i — I — i — 



50 100 150 200 250 300 350 



Figure 1: Numerical solutions of Eq. (13. 7p obtained by using the three maps. Note 
that b = 90° — i?, and I = ip, where (b, I) are the galactic coordinates. 



In Fig. 1, we plot the numerical solutions d versus p (which define the direction 
of the symmetry axis), using the so-called galactic coordinates system characterized 
by the galactic latitude b, and galactic longitude I. In our notation, the angle b 
corresponds to b = 90° — $ while / = p. From Fig. 1 we see that the galactic latitude 



11 



of the symmetry axis is remarkably independent on the adopted CMB temperature 
fluctuation map. On the contrary, we find that the galactic longitude is poorly 
constrained, for it may vary in a large interval. Indeed, we have b ~ 50° — 54°, while 
40° < / < 140°, 240° < / < 310°. One could think of choosing the value of<p (for each 
map) by imposing that the corresponding solution (ed ec , <p, 02, 03) has the smallest 
difference between (02, 03) and those in absence of anisotropic component listed in 
Table 1. Since such a difference indicates how much the intrinsic quadrupole needs 
to be rotated to accommodate the anisotropic component, the above choice would 
represent the "least radical modification" to the standard (isotropic) cosmological 
model. For the three maps SILC400, WILC3YR, TCM3YR, we could then pick the 
solutions corresponding to I — <p ~ 255°, 88°, 83°, respectively. Also in this case, 
however, we cannot univocally fix the galactic longitude of the magnetic field. 

It is interesting to stress that our determination of the direction of the symmetry 
axis is in fair agreement with the statistical analysis of Ref. [21] which confirmed 
the strong alignment and planarity of the quadrupole and the octupole. Thus, we 
see that our proposal for an ellipsoidal universe could provide a natural solution to 
the low-£ CMB puzzles. 

4 Limits on Cosmic Magnetic Fields from Big Bang Nucleosynthesis 
and Large Scale Structure formation. 

In the previous Section, we have seen that ellipsoidal universe with an eccentric- 
ity at decoupling ed ec — 0.64 x 10~ 2 could provide a natural solution to the low 
quadrupole CMB puzzle. We believe that the most interesting and intriguing pos- 
sibility is plane-symmetric geometry induced by cosmological magnetic fields. In 
fact, magnetic fields are observed in the universe up to cosmological scales. So that, 
we cannot exclude that the origin of the presently-observed cosmic magnetic fields 
is primordial as long the predictions of the standard cosmological model are not 
invalidate. Indeed, our estimate for the present-time magnetic field strength, 



is in agreement with the observed magnetic fields within galaxies and clusters of 
galaxy [13]. However, the presence of a cosmic primordial magnetic field must also 
fulfil the constraints coming from Big Bang Nucleosynthesis (BBN) and Large Scale 
Structure (LSS) formation. 

Limits from BBN. Since in the early universe the conductivity of the primordial 
plasma is very high, magnetic fields are frozen into the plasma and evolve adia- 
batically, B oc a -2 , where a oc g^^T^ 1 , T being the temperature, and g*s counts 
the total number of effectively massless degrees of freedom referring to the entropy 
density of the universe [26J. 

The limit coming from BBN refers to uniform magnetic fields at that time. The 
upper bound is given in Ref. [2"T] : 



B 



4.6 x KT 9 G 



(4.1) 



B(T BBN ) < 1 x 10 n G 



(4.2) 
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where T r 



BBN — 



10 9 K 



O.lMeV. This limit translates into: 




2/3 



) 



2 



B 



o — 



B(T BBN ) < 6 x 1(T 7 G, 



(4.3) 



BBN 



where B = B(T ), and we used c/* 5 (T BBN ) ~ g*s(T ) ^ 3.91, and T ~ 2.35 x l(T 4 eV 
[26J. We see, indeed, that our estimate (14.11) does not violate the upper bound (14. 31) . 

Limits from LSS. Anisotropic cosmological models have been extensively studied 
since long time (see, for instance, Ref. [28] and references therein). It is known that 
strong anisotropic expansion could influence the formation of large-scale structures. 
In particular, in case of strong anisotropies in the expansion of the universe the time 
of growth of structures could increases by a factor 3 — 5 [28J . 

In matter-dominated universe, we see from Eq. (12.191) that the eccentricity evolves 
as e 2 (z) oc z 3 ^ 2 . So that, at the epoch of matter-radiation equality we get: 



where e eq = e(z eq ) and we used Zdec — 1088, z eq ~ 3400 [26], and ed ec — 0.64 x 10~ 2 . 

Thus, we see that at the epoch of matter-radiation equality the anisotropy in the 
cosmological expansion is small. So that, we do not expect dramatic effects on 
the processes of growth of structures. Nevertheless, it could well be that small 
anisotropies in the cosmic expansion could leave signatures on large scales which 
should be contrasted with observations. 

5 Conclusions 

The recent measurements of the cosmic microwave background angular spectrum 
have greatly improved our understanding of the universe, in particular we have 
information on the properties and origin of the density fluctuation of the cosmic 
plasma. In fact, the observed CMB fluctuations are in remarkable agreement with 
the prediction of the ACDM standard model with scale-invariant adiabatic fluctua- 
tions generated during the inflationary epoch. However, some anomalies have been 
discovered on large angular scales including, in particular, the low CMB quadrupole 
power and the alignment and planarity of the quadrupole and octupole modes. It 
should be stressed, however, that the large-angle anomalies in the CMB anisotropy 
are still subject to an intense debate [29J. For instance, it has been suggested that 
the CMB anomalies could be explained by the residual Galactic foreground emis- 
sion [30]. Therefore, it is very important to improve and develop techniques which 
allow a better removal of residual foreground contamination [31 J . 

In this paper, we have further elaborated our previous proposal [10] that an "el- 
lipsoidal expansion" of the universe could resolve the CMB quadrupole anomaly. We 
have shown that such anisotropic expansion (described by a plane-symmetric met- 
ric) can be generated by cosmological magnetic fields or topological defects, such as 
cosmic domain walls or cosmic strings. Indeed, topological cosmic defects are relic 




(4.4) 




13 



structures that are predicted to be produced in the course of symmetry breaking in 
the hot, early universe. Nevertheless, we believe that the most interesting and in- 
triguing possibility is plane-symmetric geometry induced by cosmological magnetic 
fields. In fact, magnetic fields have been already observed in the universe up to 
cosmological scales. 

We have shown that the quadrupole anomaly can be resolved if the last scat- 
tering surface of CMB is an ellipsoid. Indeed, we found that, if the eccentricity at 
decoupling is: 

e dcc ~ 0.64 x 10" 2 , (5.1) 

then the quadrupole amplitude can be drastically reduced without affecting higher 
multipoles of the angular power spectrum of the temperature anisotropy. Remark- 
ably, our estimate of ed ec gives for the strength of the cosmic magnetic field: 

B ~ 4.6 x 1(T 9 G, (5.2) 

which agrees with the limits arising from primordial nucleosynthesis and large scale 
structure formation. Moreover, we have obtained tight constraints on the direction 
(b, I) of the axis of symmetry: 

b ~ 50° - 54°, 40° < I < 140°, 240° < I < 310°, (5.3) 

where b and / are the galactic latitude and the galactic longitude, respectively. This 
constraints are in fair agreement with recent statistical analyses of the cleaned CMB 
temperature fluctuation maps of the WMAP3 data obtained by using an improved 
internal linear combination method as Galactic foreground subtraction technique. 

In conclusion, our proposal for the ellipsoidal universe offers a natural solution to 
the CMB quadrupole anomaly which future pattern searches with more refined data, 
such as further WMAP data releases or PLANK data, will be able to confirm or 
reject. In addition, recently it has been shown [32] that the large scale polarization 
of the cosmic microwave background induced by an ellipsoidal universe compares 
quite well to the average level of polarization detected by the Wilkinson Microwave 
Anisotropy Probe. Still, it is necessary to better understand the foreground con- 
tamination of the polarization to reach a firm conclusion. 

Finally, we find amusing that there are already independent indications of a sym- 
metry axis in the large-scale geometry of the universe, coming from the analysis of 
spiral galaxies in the Sloan Digital Sky Survey [33] and the analysis of polarization 
of electromagnetic radiation propagating over cosmological distances [31] . 
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6 Appendix 

In this Appendix, we solve the system of equations (13 .7p for 1 = 2. The numerical 
values of the a 2m 's are listed in Table 1, while the a^ m 's and aJ, m 's are given by 
Eq. (13.141) and Eq. (13. 16j) . respectively. 

Since it is always possible to choose 020 real, we take 0i = 0. Moreover, the 
temperature anisotropy is real so that we have a 2 - m = (— l) m (a 2m )*, a 2- m = 
(— l) m (a^ m )*, and a\_ m = (— l)' m (a 2m )*. Therefore, the system of equations (13. Tj) 
reduces to: 

a 20 = --^=[l + 3cos(2tf)]eL + Q h (6.1) 



6v/5 V n dec V 3 



/ 7T / 7T 

Re[a 2 i] = J — cosip sin(2$) e 2 dcc + J - cos0 2 Qi, (6.2) 

/ 7T / 7T 

Im[a 2 i] = - W — siny? sin(2$) e 2 dcc + J- sin0 2 Qi, (6.3) 

71 ■ 2„a „2 ^ 



Re[a 22 ] = - J — cos(p sin $e dcc + J - cos0 3 Qi, (6.4) 

Im[a 22 ] = W— sinyj sin^e 2 ^ + J - sin0 3 Qi, (6.5) 

where Qi is given by Eq. (13.171) . We see that these equations form a system of 5 
transcendental equations containing 5 unknown parameters: ed ec , f, 02, and 03. 
Solving Eq. (16.11) with respect to $ we get two independent solutions: 

tf = {tf,7T -tf}, (6.6) 

where 



~ 1 ( 5y/5nQi - 5\/l5a 20 - V^tt e 2 Acc \ 

^=-arccos . (6.7) 

2 V 3 V37re2 ec y 

Squaring Eqs. (16.41) and (16.51) . adding side by side, and then solving with respect to 
(p, we obtain 8 independent solutions: 



where 



and 



^ = 2 arCCOS I > (6-9) 



2ir 

a = \ — Re[a 22 ] sin 2 ?9e 2 lec , (6.10) 

V -LO 

(3 = ^/-Im[a 22 ]sin¥eL, (6.11) 

7 = ^sin^eL + M 2 -^: 2 . (6.12) 
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By dividing side by side Eqs. (16.31) and (16. 2p . and solving with respect to 2 , we get: 



tan02 = v 30Ma 21 ] + vising sin(20) ^ 

dec 



30Re[a 2 i] — y^sin^ sin(2-#) e 
The same procedure applied to Eqs. (16.51) and f !6. 31) results in: 

= v^ImM-Vi^) sin^eL 
v^0Re[a 22 ] + 0Fsim» sin 2 tfe 2 ec 

Finally, by squaring Eqs. (16.21) and (16. 3p . and adding side by side, we get: 



e% x +2cei x -d = 0, (6.15) 



where we have defined 



/30 

c((p,$) = \ — (Re[a 21 ] cos^ — Im[a 2 i] siny?) csc(2^), (6.16) 

V 7T 

d((p,&) = lO^Q 2 - ^ |a 21 | 2 ^csc 2 (2tf) . (6.17) 

We observe that the couple (■#, <p) can assume 16 different values, according to 
Eqs. (16.61) and (16. 8p . Inserting these values in Eqs. fl6.15[) -( T67TT|) we arrive at 16 
different equations for ed ec - It is straightforward to verify that only 8 of these are 
"independent", in the sense that, given a solution (ed ec , V 9 ) °f one °f t ne "indepen- 
dent" equations, then (edec, 7 !" — V 9 =t 7r ) is a solutions of one of the "dependent" 
ones (we must take the plus sign if ip < it and the minus sign if ip > 7r). The 8 
independent equations can be solved numerically and their solutions are presented 
in Table 2, 3, and 4. Here, we just observe that, to the map SILC400 it corresponds 
8 independent solutions, while to the maps WILC3YR and TCM3YR there corre- 
spond only 4 independent solutions. 

Finally, we derive the approximate solutions (I3.20p - (l3.2ip . To this end, we may 
formally solve Eq. (I6.15P to get e 2 ec = c ± Vc 2 + d. If Q\ ^> |a 2 i| then d ^> c 2 , 
and we obtain e\ ec ~ Vl~0 Qi | csc(2$)|. Now, it is easy to check numerically that 
|csc(2#)| ~ 1. Indeed, from Table 2, 3, 4, we get 1.01 < |csc(2#)| < 1.07. As a 
consequence, we have e 2 dec ~ vlOQi, which indeed agrees with Eq. (13.201) . Inserting 
Eq. (|MDD into Eq. (JEU) we easily recover Eq. (EEZIjl . 
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